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Diffraction in time of a particle confined in a box which its walls are removed suddenly at t = 
is studied. The solution of the time-dependent Schrodinger equation is discussed analytically and 
numerically for various initial wavefunctions. In each case Bohmian trajectories of the particles are 
computed and also the mean arrival time at a given location is studied as a function of the initial 
state. 
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I. INTRODUCTION 

Several works have been done on time-dependent boundary conditions [H-@|. Diffraction in time was initially 
introduced by Moshinsky [l[ . A beam of particles impinging from the left on a totally absorbing shutter located at 
the origin which is suddenly turned off at time zero. The transient current has a close mathematical resemblance 
with the intensity of light in the Fresnel diffraction by a straight edge. An interesting feature of the solutions for 
cutoff initial waves, occurring both in the free case [l8| and in the presence of a potential interaction is that, if 
initially there is a zero probability for the particle to be at x > 0, as soon as t = + , there is instantaneously, a finite, 
though very small, probability to find the particle at any point x > 0. This non-local behavior of the Schrdinger 
solution is due to its nonrelativistic nature and not a result of the quantum shutter setup Application of the 
Klein-Gordon equation to the shutter problem [l[ shows that the probability density is restricted to the accessible 
region x < ct (c is the speed of light). See [l(| for a recent review. Gerasimov and Kazarnovskii confined the 
initial wave in a finite region by introducing a second shutter at the point x = L. Godoy (ri| pointed out the analogy 
with Fraunhoffcr diffraction in the case of small box (compared to the de Broglie length), and Fresnel diffraction, 
for larger confinements. In this context, by considering the problem of a particle in a one-dimensional box potential 
which its walls are suddenly removed at some time, the aim of the present paper is to probe some aspects of the 
time-dependent boundary condition for a particle confined in an square well focusing on Bohmian interpretation of 
quantum mechanics that have remained hitherto unnoticed. The computed Bohmian trajectories are instructive in 
revealing the conceptual ramifications of such an example. 

Although the formalism of Bohmian mechanics does not give predictions going beyond those of QM whenever the 
predictions of the later are unambiguous, it should be favored because of its interpretational advantages stemming 
from the ontological continuity between the classical and the quantum domains [131 ]. Noticing the Bohmian arrival 
time formulation by means of cut-off current, it has been argued predictions of Bohmain mechanics are in contradiction 
to the standard quantum mechanical formalism fl4| . In nonrelativistic Bohmian mechanics the world is described by 
point-like particles which follow trajectories determined by a law of motion. The evolution of the positions of these 



particles are guided by a wavefunction which itself evolves according to the Schrodinger equation [15l4l9| . In this 
5^ ■ theory, in the absence of any measuring device, one finds [20r]22j that for those particles that actually reach x = X, 
the arrival-time distribution is given by the modulus of the probability current density, i.e., \j(X, t)\. We will proceed 
as follows. In Sec. HI] solution of time-dependent Schrodinger equation is given for a particle which is initially confined 
in a box. Sec. IIIII contains a very brief review of relevant parts of Bohm's interpretation of quantum mechanics. Sec. 
IIVI gives numerical results. Finally, in Sec. |V]wc present the concluding remarks. 
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II. FREE PROPAGATION OF A PARTICLE INITIALLY CONFINED IN AN SQUARE WELL 

Consider a particle which is initially confined in an interval [0,L] with wavefunction tpo(x). If at time t = it is 
free, then at any instant t > its wavefunction is given by, 



ip(x,t) = / G(x,t\x' ,0)ip (x') dx', (1) 

J — OO 

in which G(x, 0) is the free particle propagator and is determined by, 



G(x,t\x',0) = ^ei^-'> , (2) 

and ipo i x ') is t ne initial wavefunction. In this work we take initial wavefunction to be, (a) a stationary state of a particle 
inside a well with hard (perfect reflective) walls at x = and x = L and (b) a motionless localized Gaussian wave- 
packet in that region with negligible overlap with the walls of the well. To avoid any problem concerning the boundary 
conditions, one can suppose in this case that the walls act as absorbers or the tails of the wave packet has been cut by 
the walls of the well. In the first case initial wavefunction is given by, ipo{x) — (f> n (x) = y/2/ L sm(k n x) X[o,l](%), with 
k n = nir/L. X[o,L]( x ) — &(L — x) — &(—x) is the characteristic function in the interval [0,L]. Simultaneous removal 
of both walls leads at time t to the wavefunction [Tl], Qjl, HH 



L 



ip n (x,t) = a/^^YZ I X)2 sin (k n x') dx' , 



= 1p n!+ (x,t) + 1p n -(x,t) , (3) 

which is a superposition of a right and a left moving diffracted in time plane waves. After doing some simple algebra, 
one gets 

= -^=e ik ^- iE ^ h [F n {x -L,t)- F n (x,t)} , (4) 

TP n ,-(x,t) = -jl=e- ik " x - iE ^ h [F n (L-x,t)-F n (-x,t)] , (5) 
V4i 3 L 



with E n = h 2 k^ l /2m and 



F n (x,t) = / due 1 ™ I 2 , (6) 



with upper limit £ n (x,t) = ^jjj (v n t — x) in which v n = hk n /m. Let us for later use, compute the derivative of 
ip n {x,t) with respect to x. 

gj/VhOM) , , ... ik n x-iE n t/n ^F n {x - L,t) dF n (x,t) 

^- — iK n ip„ t+ {x, i) -+- e [ — — J , ( I) 

g%M = ^^_ {Xit) + e ^^t/H [ dFn(L-x 1 t) _ dF n (-x,t) 
ox ox Ox 

in which, 



dx V 7T?lt 

Now, by some straightforward algebra, one can show, 



dF n (x,t) (g) 



^L =L/2 = 2e'^ h co S (k n L/2) ([F n {-L/2,t) - F n (L/2,t)] + ^[e^^ 1 -e ^^ ^j , (10) 
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which is zero for odd n. Note that one can find this, without doing any algebra. Wavcfunction is an even (odd) 
function for odd (even) n with respect to the point x = L/2, so its derivative is an odd (even) function for odd (even) 
n with respect to that point. 



In the second case ipo(x) = j2~2yj& e 4 "° X[o,l](%), hi which xq is the center of the packet and cto is its rms 
width; cr = (x 2 ) — (x)q. After simultaneous removal of both walls, wavcfunction is given by, 



I m f L , - <x ~\ 0) +igi.( x - x >) 2 



III. BOHMIAN TRAJECTORIES 



In nonrelativistic Bohmian mechanics the world is described by point-like particles which follow trajectories de- 
termined by a law of motion. The evolution of the positions of these particles are guided by a wavcfunction which 
itself evolves according to the Schrodingcr equation. Given the initial position x' ' = x(t = 0) of a particle with 
the initial wavcfunction tpo(x), its subsequent trajectory x{x^\t) is uniquely determined by simultaneous integration 
of the time dependent Schrodinger equation, and the guidance equation dx j^ = v(x(i),t), in which v = i, where 

j = -j^Q 1S the probability current density and p = \tp(x,t)\ 2 is the probability density. In the context of 

Bohmian mechanics arrival time distribution at a given location, say x = X is given by f20l - |22j . 



Jo dt\j{X,t)\ 

So, mean arrival time at observation point x — X is determined by, 

POO 

t(x) = dtt n x (t) . (13) 

Jo 

A general formulation for Bohmian arrival times was given in [23| and a formula for the numerical calculation of such 
Bohmian arrival times in the case of ID rigid inertial detectors (exactly the cases we have here) was presented in [24| . 
The derived formula there doesnt require the explicit calculation of the Bohmian trajectories and the resulting 'cut-off 
current' can be considered to be a generalization of the arrival time probability density introduced by Leavens (20l - [2^ [. 



According to Eq. (12) of [24j . the probability density ILx(t) of the arrival time distribution for a point detector at 
x = X takes the form, 

H x (r) - (MP(t)) 1 j(A,r)[e(/x(r)-max{/ x ( S )/0< S <r})-e(-/x(T)-max{-/ x (.s)/0< S <r}X]l4) 
in which P is the detection probability, 

P(t) = max{/ x (s)/0 < s < t} + max{-/ x (s)/0 < s < t} , (15) 

with 

f x (s) = [ j(X,t)dt. (16) 



For the case of positive or negative j(X,t) Eq. (fT4)) reduces to Eq. (fT2|) . 



IV. NUMERICAL RESULTS 



For numerical calculation width of the well is chosen as L = 1 /im. All of the calculations are presented for 
Rubidium atoms with mass m = 1.42 x 10~ 25 kg. Fig. [1] shows probability density versus distance i(/im) for state 
n — 6 at different times. At longer times there are two spatial packets placed about the box and moving apart. As 
pointed out by del Campo and Muga [l2| this takes place after the semiclassical time t„ = mL 2 /2mrh as a result 
of the mapping of the underlying momentum distribution to the density profile expected asymptotically. For our 
parameters t n = (0.214/n), hence, = 0.031 ms and £500 = 4.28 x 10~ 4 ms. Fig. [5]shows probability current density 
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as a function of time at observation point x — 2 /im outside the box after removal of the walls for various stationary 
states. From Eq.([3]) it is obvious that x = L/2 remains a node of the wavefunction for even n, i.e., tp n (x = L/2, t) = 
for even n. So because of a well-known property of Bohmian paths, Bohmian particle cannot be initially located at 
x = L/2 or even pass through this point for even values of n. As mentioned above, right after Eq. (fTDj) . for odd n 
distance-derivative of the wavefunction is zero at point x = L/2. Thus, current probability density and consequently 
Bohmian velocity is zero in this point all the time. Therefore, a Bohmian particle which is initially at x^ = L/2, 
will remain at rest. Because of noncrossing property of Bohmian path, particles with a;(°) < L/2 (x^ > L/2) will 
go backwards (forwards). Sec Fig. [31 We have used Runge-Kutta method for the simultaneous integration of the 
time-dependent Schodingcr equation and the guidance law to compute Bohmian trajectories. For n = 7, after t = 0.03 
ms trajectories exhibit a bifurcation into two main branches while for n = 500 this takes place after t = 0.0004 ms. 
These values coincide very well with tj and t 5m in the above. Fig. U shows the mean arrival time, in the context 
of Bohmian mechanics, as a function of quantum number n at the observation point x = 2 (jxn. It is clear that r 
decreases with n as one expects, because by growing n, semiclassical velocity increases. 

In the case of a Gaussian packet, parameters of the packet are chosen as x$ = 0.5 /j,m and ao = 0.25 ^m. It should 
be noted that with these parameters, initial Gaussian packet is not normalized to unity but to 0.954543 (truncated 
Gaussian packet). If some one choses initial packet narrower than us, in such a way that it locates totally inside the 
well, then after removing the walls its evolution will be the same as a free Gaussian packet which is not desired here. To 
show the differences we consider a free Gaussian packet with the same parameters as well. Fig. [5]shows the probability 
density versus distance, 0.1 ms after removal of the walls and the probability current density at observation point 
x = 2 /im as a function of time for a free and a confined truncated Gaussian packet. In the confined case one sees some 
oscillations in the plot of current density which arc absent in the free case. Finally, Fig[B]shows a selection of Bohmian 
paths for both cases. In the free case trajectories are determined by x(t) = xq + (x^ — xo)yl + (Tit / '2ma^) 2 , where 
is the initial position of the particle [l8j. From the figure it follows that the trajectory which starts at x^ = x is 
the bifurcation trajectory in both cases and Bohmian velocity of a path in confined case is larger than the Bohmian 
velocity of the corresponding path in the free case. 

V. SUMMARY AND CONCLUSION 

The dynamics of particles, with various initial wavefunction, released from a box has been studied. Different 
boundary conditions, the absorbing wall and the reflecting one, is considered. Such studies suggest that the time- 
varying boundary conditions can give rise to interesting action-at-a-distancc effects in quantum mechanics. Quantum 
temporal oscillations of matter waves released from a confinement region constitute the hallmark of the diffraction in 
time effect. Mean arrival time at a observation point outside the box has been considered for various initial states. 
Moreover, Bohmian paths of the particle are computed. Our calculated mean arrival time may have no relevance 
to the experiment. Within Bohm's causal theory of quantum mechanics if we try to measure properties other than 
position (the only intrinsic property), we find that the result is affected by the process of interaction in a way that 
depends, not only on the total wave function, but also on the details of the initial conditions of both the particle and 
the apparatus. Our calculations can be verified if mean arrival times can be experimentally measured in such a way 
that the experiment does not perturb the unmeasured quantity. 
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FIG. 2: Probability current density (1/ms) as a function of time t (ms) at observation point x = 2 /im for states (a) n = 1, (b) 
n = 50, (c) n = 100 and (d) n = 150. 




FIG. 4: Mean arrival time at detector position x = 2 /xm for different states. 
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x (Jim) t (ms) 

FIG. 5: Probability density (1/^m) versus distance x(fim) at time t = 0.1 ms for (a) a free motionless Gaussian wave-packet 
and (c) a motionless truncated Gaussian wave-packet initially confined in a box. Probability current density (1/ms) versus time 
t (ms) at observation point x = 2 /im for (b) a free motionless Gaussian wave-packet and (d) a motionless truncated Gaussian 
wave-packet initially confined in a box. 




FIG. 6: A selection of Bohmian paths for (a) a free motionless Gaussian wave-packet and (b) a motionless truncated Gaussian 
wave-packet initially confined in a box. 



